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ABSTRACT 

Context. CoRoT-2b is one of the most anomalously large exoplanet known. Given its high mass, its large radius cannot be explained 
by standard evolution models. Interestingly, the planet's parent star is an active, rapidly rotating solar-like star with with spots covering 
a large fraction (7-20%) of its visible surface. 

Aims. We attempt to constrain the properties of the star-planet system and understand whether the planet's inferred large size may be 
caused a systematic error in the inferred parameters, and if not, how it can be explained. 

Methods. We combine stellar and planetary evolution codes based on all available spectroscopic and photometric data to obtain self- 
consistent constraints on the system parameters. 

Results. We find no systematic error in the stellar modeling (including spots and stellar activity) that would cause a ~ 10% reduction 
in size of the star and thus the planet. Two classes of solutions are found: the usual main-sequence solution for the star yields for 
the planet a mass of 3.67 + 0.13 Mj, a radius of 1.55 + 0.03 Rj for an age that is at least 130 Ma and should be younger than 500 Ma 
given the star's rapid rotation and significant activity. We identify another class of solutions on the pre-main sequence, for which the 
planet's mass is 3.45 + 0.27 Mj and its radius is 1.50 + 0.06 Rj for an age of 30 to 40 Ma. These extremely young solutions provide 
the simplest explanation of the planet's size that can then be matched by a simple contraction from an initially hot, expanded state, 
if the atmospheric opacities are larger by a factor of ~ 3 than usually assumed for solar composition atmospheres. Other solutions 
imply that the present inflated radius of CoRoT-2b is transient and the result of an event that occurred less than 20 Ma ago, i.e., a 
giant impact with another Jupiter-mass planet, or interactions with another object in the system that caused a significant rise in the 
eccentricity followed by the rapid circularization of its orbit. 

Conclusions. Additional observations of CoRoT-2 that could help us to understand this system include searches for an infrared 
excess, a debris disk, and additional companions. The determination of a complete infrared lightcurve including both the primary and 
secondary transits would also be extremely valuable to constrain the planet's atmospheric properties and determine the planet- to- star 
radius ratio in a manner less vulnerable to systematic errors caused by stellar activity. 

Key words. Star: individual: CoRoT-2; (Stars:) planetary systems; Stars: pre-main sequence; Planets and satellites: physical evolution 



1. Introduction 

CoRoT-2b is the second transiting planet discovered by the space 
mission CoRoT (Alonso et al. 2008; Bouchy et al. 2008). It is 
noteworthy because of its relatively large mass and large size. 
As such, it belongs to the class of anomalously large extrasolar 
giant planets, i.e. planets that are larger than predicted by stan- 
dard theoretical models for the evolution of an irradiated, solar- 
composition gas giant (Guillot et al. 2006). 

Anomalously large planets require additional heat sources 
(e.g. Bodenheimer et al. 2001; Guillot & Showman 2002; 
Baraffe et al. 2003; Guillot et al. 2006; Burrows et al. 2007) to 
explain their radii. It is now clear that these objects are common, 
representing about half of the presently known transiting plan- 
ets (Guillot 2008; Miller et al. 2009). However the properties 
of CoRoT-2b are probably the most difficult to explain: it has 
one of the largest radius anomalies (defined as the difference be- 
tween observed and modeled for a solar-composition irradiated 
planet), and yet is massive, meaning that modifying its evolution 
to slow its contraction is difficult and requires significant devia- 
tions from the standard models. 
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Secondary transits of CoRoT-2b have been detected by sev- 
eral instruments providing puzzling results: Measured brightness 
temperatures have been found to vary from 1325 ± 180 K at 8 fim 
to 1805 ± 70 K at 4.5 jim (Gillon et al. 2010; Snellen et al. 2010), 
and up to 2170 ± 50 K in the visible (Snellen et al. 2010; Alonso 
et al. 2010b), significantly higher than the zero-albedo equilib- 
rium temperature of the planet 1530 ± 140 K. 

The star CoRoT-2 is also remarkable because of its un- 
usual variability and spot coverage: analyses of the CoRoT 
lightcurve for this system -spanning 152 days of nearly- 
continuous observations- led to a precise determination of the 
star's rotation period, 4.52 ± 0.14 days (Lanza et al. 2009), and 
estimates of the spot coverage that range between 7% and 20%, 
for a spot contrast of between ~ 0.3 and 0.7 (Lanza et al. 2009; 
Huber et al. 2010; Silva-Valio et al. 2010). The spot cover- 
age (i.e. the fraction of the stellar surface that is occupied by 
starspots) may even locally reach 37% for a contrast of 0.7 (sim- 
ilar to the average value for sunspots) at the latitudes where the 
planetary eclipses occur (Huber et al. 2010). 

Could the peculiarities of the star account for the unusually 
large inferred size of the planet or do we require additional heat 
sources in the planet? We address this question by first revisit- 
ing the star's evolution by accounting for the presence of spots 
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Table 1. Observational constraints on the stellar parameters 





Spectroscopy 


7eff,Hor 


5450+ 120 K 


Bouchy et al. (2008) 


T c S,\Fc] 


5625 + 120 K 


Alonso et al. (2008) 




5608 + 37 K 


Ammler-von Eiff et al. (2009) 


logg 


4.3 + 0.2 cms -2 


Alonso et al. (2008) 


logg 


4.71 +0.2 cms" 2 


Ammler-von Eiff et al. (2009) 


[M/H] 


0.0 + 0.1 


Alonso et al. (2008) 




Photometry and RV 


My 


12.57 


Alonso et al. (2008) 


p+ 


4.522 + 0.024 days 


Lanza et al. (2009) 


Povb 


1.7429935 + 10" 6 days 


Gillon et al. (2010) 


k 


0.1711 + 0.0011 


This work 


k 


0.1658 + 0.0004 


Gillon et al. (2010) 


k 


0.1720 + 0.0010 


Czesla et al. (2009) 


k 


0.1667 + 0.0006 


Alonso et al. (2008) 


i 


88.08^ deg 


Gillon et al. (2010) 


P* 


1.870 + 0.026 gem" 3 


Alonso et al. (2008) 


P* 


1.814^-050 gcm -3 


Gillon etal. (2010) 


K 


603 + 18 ms" 1 


Gillon et al. (2010) 


a/R+ 


6.64 + 0.03 


Gillon et al. (2010) 



(Sect. 2), then applying these results to planetary evolution mod- 
els (Sect. 3). 

2. The evolution of a spotted star: Constraints on 
CoRoT-2 

2.1. Constraints derived from spectroscopic and transit 
photometry analyses 

The star CoRoT-2 has been identified as a G7 dwarf of solar 
composition by spectroscopic analyses. Table 1 lists its known 
physical properties inferred from spectroscopy, radial velocime- 
try, and the analysis of its transit lightcurve. 

The effective temperature of CoRoT-2 differs slightly be- 
tween measurements: Bouchy et al. (2008) report a low value 
from H a bands and relatively low signal-to-noise HARPS spec- 
tra, and a high value, also from HARPS but using Fe spectral 
lines. Ammler-von Eiff et al. (2009) essentially confirm this last 
value with Fe lines and UVES, but with a smaller error bar 
(which does not however include systematic effects). 

The spectral determination of the star's gravity is, as com- 
monly found for stars, quite uncertain. Measurements inferred 
from HARPS and UVES data differ slightly, there being one cr 
error bars that barely overlap. 

As is typical of a star with a transiting planet, the most 
stringent constraint is that on the stellar density. For CoRoT-2, 
Alonso et al. (2008) are able to determine the duration of the 
transit so precisely that the stellar density is constrained to within 
only 1.4%. 

Given these measurements, Alonso et al. (2008) infer plane- 
tary parameters, M p = 3.31 ± 0.16 Mj and R p = 1.465 ± 0.029 
Rj. This implies that CoRoT-2b is extremely inflated even rela- 
tive to other large transiting planets. Before we attempt to model 
it, we estimate the accuracy to which the stellar parameters can 
really be derived. For example, this initial estimate by Alonso 
et al. assumes a circular orbit and no influence by spots on the 
photometric determination. 

A subsequent analysis of the CoRoT-2 lightcurves (Czesla 
et al. 2009) shows that the presence of spots during the transits 
affects the photometry and in particular the depth k 2 of the tran- 
sits: when the planet transits, it occasionally occults star spots, 



thereby blocking a smaller fraction of the stellar light. On aver- 
age, the transits appear less deep than for a star without spots, 
implying that the planetary radius is underestimated when ne- 
glecting the effect of spots. 

The effect of a non-zero eccentricity is to modify the rela- 
tion between transit duration and stellar density in a non-trivial 
way relative to a circular orbit (e.g. Tingley & Sackett 2005). 
Gillon et al. (2010) refined the analysis of Alonso et al. us- 
ing constraints on the eccentricity and argument of the perias- 
tron from the radial velocity data within a Markov Chain Monte 
Carlo approach. They found that solutions with a slight eccen- 
tricity (e ~ 0.015) are more likely and thus inferred a lower 
stellar density p* and a slightly smaller value of k. 

Most of the parameters used for this work are based on the 
analysis of Gillon et al. (2010) (see Table 1). However, to ac- 
count for the effect of spots on the transit depths, and because 
the Czesla et al. (2009) study did not allow for the possibility 
of an eccentric planet, we choose to derive a probable value 
of k that accounts for both spots and a non-circular orbit by 
a simple proportionality relation between the different studies: 
k = ^czesla /^Alonso x feiiion- The error bar in k is calculated to be 
the quadratic mean between the values of k found by the Czesla 
et al. and Gillon et al. studies. 



2.2. Evolution models 

Stellar evolution models are needed to derive the star's and 
planet's masses and radii. In most of this work, we use a grid 
of quasi-static evolutions for stars with masses between 0.6 and 
1.3 M (AM* = 0.005M o ) calculated with the CESAM evo- 
lution code (Morel & Lebreton 2008). The grid has been cal- 
ibrated with respect to the Sun, which is most accurately de- 
scribed (< 10" 4 relative precision) by mass fractions of hydro- 
gen X Q = 0.7065, helium Y Q - 0.2740, and heavy elements 
Z Q = 0.0195, based on the actual Z/X = 0.0245 ratio of Grevesse 
& Noels (1993), and spectroscopic parameters T e ff, = 5778 K, 
and L = 3.846 x 10 26 W for an age of 4.57 Ga. A standard 
mixing-length approach without overshooting is used in the en- 
ergy transport equation. Our calibrated Sun has a mixing length 
parameter a Q = 2.052. The atmospheric boundary is calculated 
using a T(r) relation derived from MARCS models (Gustafsson 
et al. 2008), and we consider the microscopic diffusion of chem- 
ical elements in the radiative zone (therefore the atmospheric 
metallicity varies as a function of time). We chose the abun- 
dances of Grevesse & Noels (1993) because the seismological 
constraints are not met when using other, more recent abun- 
dances. 

For comparison, we also use the grid of models calculated by 
Baraffe et al. (1998) (hereafter BCAH98) for solar-composition 
stars. Those models assume a non-grey atmosphere, and the con- 
vection is treated in the context of the mixing-length theory with 
no core-overshooting and no diffusion. The tables assume that 
X = 0.716, Y = 0.282, Z = 0.020, and a = 1.9. For a so- 
lar model (1 M Q , 1 R Q ), the tables yield an effective temperature 
T eff ,o = 5797 K and luminosity L Q = 3.801 x 10 26 W for an 
assumed age of 4.61 Ga. 

Finally, results are also compared to the so-called Y 2 stel- 
lar evolution tracks (hereafter YY) for solar composition stars 
(Demarque et al. 2004). These models include convection- 
overshoot and diffusion. 
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2.3. The effect of spots and activity 

By definition, the effective temperature of a star is linked to 
its total luminosity L and radius R by the well-known relation 



4/r/?Vr e 4 ff 



(1) 



However, in the presence of spots and activity, this relation has 
to be revised, because the correspondence between the effective 
temperature derived from a spectroscopic analysis (which we re- 
fer to as r e ff ? spectro) and that obtained from Eq. 1 no longer holds. 

Across the face of the Sun, spots are indicative of high mag- 
netic activity: although the spots block a fraction (~ 0.1%) of 
the starlight, their appearance is connected to a global increase 
in the total solar irradiance by about 0.1%, from ~ 1365.5 to 
1366.5 W m~ 2 due to the presence of bright faculae (e.g. Frohlich 
& Lean 2004; Krivova & Solanki 2008). This implies that the 
flux emitted in the visible and most importantly in the UV in- 
creases. Secular variations in the total solar irradiance based on 
models and a ~ 300 year record of the solar activity also amount 
to about 0.1% (Solanki & Fligge 2000). 

For stars with activity levels comparable to that of the Sun or 
lower, this level of uncertainty is much smaller than that attain- 
able from spectroscopic measurements. In the case of CoRoT-2, 
the area covered by starspots is 70 to 200 times larger than for the 
active Sun and thus can potentially affect what may be inferred 
about the star properties, i.e. its luminosity, radius and mass. 

In Table 1, the effective temperature of CoRoT-2 was de- 
termined from either the H-a line at 656.3 nm or fits in the 
visible (450 to 740 nm). For comparison, between minimum 
and maximum of activity, the Sun increases its relative flux 
AF A /F A by about 0.1%, 0.05%, and 0.07% at wavelengths in the 
range of 400-500 nm, 500-700 nm, and 700-800 nm, respectively 
(Krivova et al. 2006). Variations in the brightness temperatures 
at these wavelengths (proportional to crT 4 ) are 4 times smaller. If 
we assume that these number indeed scale with the surface area 
covered by spots (a big if arguably), the potential mismatch be- 
tween the spectroscopically inferred r e ff ? spe ctro and the true 
amounts to up to 5%. 

To see how the results may be affected by stellar activity, 
we first compare standard evolution models to evolution models 
calculated with CESAM and a modified atmospheric boundary 
condition 

L=(l-^ s )4;r/?Vr e 4 ffiSpectro , (2) 

where Xs is a factor that accounts for emission from only part of 
the surface of the star. For purely black spots with no faculae it 
corresponds to the surface fraction of spots. 

Figure 1 shows the result for a given luminosity and age, 
a star with spots simply has a higher "spectroscopic" effective 
temperature r e ff ?spectro as inferred from spectroscopic measure- 
ments than a star without spots. The evolution is simply dis- 
placed to the left of the HR diagram by an almost constant ratio 
(1 -^ s ) 1/4 in r ef f. Quantitatively, the mean deviation amounts to 
5 x 10" 5 on r ef f, with a standard deviation cr = 4.6 x 10" 4 r e ff 
and a maximum deviation 6.7 x 10" 3 r e ff. The small departures 
from this constant displacement are due to slight modifications 
of the atmospheric properties (opacities) with temperature, but 
they can safely be neglected in this study. 

When considering a star's evolution, that a star has spots 
is therefore equivalent to an added uncertainty in the mea- 
surement of its effective temperature. Activity also has the 
same consequence because it implies that both the present 
spectroscopically-determined effective temperature and the 
present luminosity may differ from their value averaged over one 




-0.5 



3.7 3.65 

log(T eff ) 



Fig. 1. Left panel: Hertzsprung-Russell evolution tracks for 
0.97 M©, Z = 0.02, a - 1.9 for a standard evolution model (red 
line), and when modifying the atmospheric boundary condition 
to account for the presence of x& = 20% of dark spots on the 
stellar surface (blue line). Right panel: differences in effective 
temperature for a given luminosity between the model with no 
spot, and the model with spots after the effective temperature has 
been shifted by a constant factor (1 -^ s ) 1/4 - 



stellar magnetic cycle. A larger error bar in T e s may be used as 
a proxy for the added uncertainty due to activity and starspots. 
Although this uncertainty may be either positive (if the luminos- 
ity of CoRoT-2 is lower than its mean value and/or the contribu- 
tion of faculae is larger than that of spots) or negative (if CoRoT- 
2 is more luminous than average and/or the spots are dark) we 
choose to only study the latter possibility. As shown in the fol- 
lowing section, only lower T e s values can yield smaller radii for 
the star and the planet and hence help to solve the inflation puz- 
zle. 



2.4. Constraints on the star's physical parameters 

The star's physical parameters (M*, R+, age) are obtained by 
matching the constraints from Table 1 to a grid of evolution mod- 
els, as depicted in Fig. 2. The two most important constraints 
of the problem are the star's effective temperature and density. 
A third constraint (not shown on the plot) is the star's present 
metallicity, which should be compared to the one obtained from 
the evolution models that include diffusion. One could include 
other constraints (such as that on log g), but in the present case, 
they are too weak to be useful. 

The quality of the fit of any given model is given by its dis- 
tance n^ to the ellipsoid of constraints, measured in units of the 
standard error in the constraints given by 



n (T = 



1/2 



(3) 



where X t are the constraints and <x; their standard deviation (as- 
sumed Gaussian). The ellipsoid of constraints (of dimension 2) 
is centered on (r ef f ,p*) and has semi-minor and semi-major axes 
km^rio-io-TtffiO-pJ, respectively, where k 2 D,n (r is the quantile of 
a 2D gaussian law at the equivalent level of confidence n a . In 
addition, &2D,iv ~ 1-52, 2.49, 3.44 for n a = 1,2,3 respectively. 
This normalization ensures that our solutions at 1,2,3 cr have 
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3.76 3.72 3.68 3.64 

log(^eff) 

Fig. 2. Evolution tracks for stars of masses between 0.85 and 
1.0M© in the effective temperature vs. stellar density space. The 
lcr observational constraints are shown as boxes: the box to 
the left (blue) corresponds to solutions that neglect the effect of 
spots. The two boxes (blue and yellow) correspond to the added 
uncertainty obtained when assuming that from to 20% of the 
stellar surface is covered by spots. CESAM tracks (plain) are 
compared to BCAH98 tracks (dashed). The models assume so- 
lar composition. Dots on the CEPAM tracks correspond to the 
following ages (small to large circles): 20, 32, 100, 500, 1000, 
5000 Ma, respectively. 

the correct probability of occurrence. However, for the metal- 
licity we adopt a relatively crude simplification: we consider as 
valid only solutions for which the [Fe/H] value is within l<T[Fe/H] 
of the measured one. We tested that in the particular case of 
CoRoT-2, considering yet larger errors in [Fe/H] has a negli- 
gible effect. In the remainder of the paper, we present models 
for which < 1,2,3, corresponding to confidence levels of 
68.3%, 95.4%, and 99.7%, respectively. 

Two possible values of the effective temperature are used: 
(i) in the no-spot case, we assume that obtained by Ammler-von 
Eiff et al. (2009) but with a slightly larger error to account for 
possible systematic errors (7^ S pectro — 

5608 ± 80 K); (ii) when 
including the effect of spots, we define a new temperature and its 
associated error A reff to take into account the presence of spots 
(from 0% to 20% of the area of the star) 

T _ T l + d-^) 1/4 

* eff — J- eff , spectro > 
Ar eff (^£r) = ^eff, spectro ~ ^eff + n^Cr r eff , spectro • 

We thus derive T eff = 5224-5688 K at lcr, T eff = 5144-5768 K 
at 2<r, and = 5064 - 5848 K at 3<x. The errors are thus in- 
trinsically non-Gaussian. 

The constraints used for the stellar density and metallicity 
are those derived by Gillon et al. (2010) and Alonso et al. (2008), 
respectively (see Table 1). 

Figure 2 shows evolution tracks in the (log T e ff, logp*) plane 
for solar-composition stars of between 0.85 and 1.0 M Q . For 
about 100 Ma, the stars contract and heat up from a low- 
temperature, low-density initial state, the so-called pre-main- 
sequence phase (PMS). They then reach a maximum in their 
density and very gradually expand while on the main- sequence 
(MS). In the case of CoRoT-2, the observational constraints on 
r e ff and p* are met either on the PMS, for ages ~ 30 Ma, or 
at much older ages in the MS phase. Only stars with masses 




3.76 3.72 3.68 3.64 

log(T eff ) 

Fig. 3. CEPAM evolution tracks for a 0.95 M solar-composition 
star for three values of the mixing length parameter around the 
Sun calibrated value (O.85a , a Q , lA5a Q ) in the effective tem- 
perature vs. stellar density space. The symbols are as in Fig. 2. 




3.76 3.72 3.68 3.64 

log(T eff ) 

Fig. 4. CEPAM evolution tracks for a 0.95 M© star for three val- 
ues of the metallicity ([Fe/H] =-0.2, 0, and 0.2, respectively) in 
the effective temperature vs. stellar density space. The mixing 
length parameter has been calibrated to the solar model. The 
symbols are as in Fig. 2.) 

between roughly 0.84 and 1.04 M intercept the box of con- 
straints at some point in their evolution. Figure 2 also shows 
that the presence of spots yields solutions at smaller masses than 
when spots are not taken into account. The agreement between 
CESAM and BCAH98 models is generally very good, with dif- 
ferences in r e ff of generally ~ 1% or less, and differences in p+ 
that can reach ~ 10% but in only a limited region of the PMS 
evolution phase. 

In Fig. 3, we explore the effect of a "reasonable" (±15%) 
modification of the mixing length parameter a on the evolu- 
tion tracks. The effect is not negligible in terms of its impact on 
both the effective temperature and the stellar density. A higher 
value of a implies a more efficient energy transport and therefore 
higher effective temperatures and generally a faster evolution. 
Less intuitively perhaps, it leads to a higher maximal stellar den- 
sity at the early stages of the MS phase. However, we emphasize 
that these models with modified mixing lengths have not been 
calibrated and do not properly reproduce the present Sun. 

The consequences of metallicity variations are shown in 
Fig. 4. An increase in the [Fe/H] value by a factor 1/3 leads to a 
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global decrease in the effective temperature by about 2%, larger 
than the lcr error in the measurements, but slightly smaller than 
the uncertainty in T q q obtained when including spots. This ef- 
fect is thus significant and shoud be included in the search for 
solutions matching the observational constraints. 

We present in Fig. 5 the ensemble of solutions in the 
age) space obtained with various assumptions. The top pan- 
els correspond to our preferred solutions using our calibrated 
CESAM evolution model and including all metallicities. At 3<x, 
a wide range of solutions is found that extends from ages be- 
tween 30 Ma and more than 10 Ga. Within l<x, two classes of 
solutions are found, either on the PMS (30 - 40 Ma) or on the 
MS (for ages > 800 Ma with no spots, or > 100 Ma when includ- 
ing spots). The range of stellar radii (which are directly propor- 
tional to the planetary radii that are inferred) extend from 0.88 
to 0.93 R at l<x, but the smallest values are obtained for either 
the youngest or oldest solutions. 

The other panels in Fig. 5 highlight the consequences of the 
different hypotheses on the solutions, when considering only so- 
lar composition models. Varying the mixing length parameter 
has consequences for the lcr solutions: a lower a value leads to 
a wider range of solutions at intermediate ages, while a higher 
a increases the separation between the very young and the very 
old solutions. However, when considering the global 3<x enve- 
lope and accounting for spots, the solutions are very similar. 

The solutions obtained by using the p* value of Alonso et al. 
(2008) are quite similar to the nominal ones, but are regarded 
as slightly over-constrained due to the assumption of a circular 
orbit. 

The last four panels in Fig. 5 provide another test of the ro- 
bustness of the solutions by a comparison with YY and BCAH98 
evolution models. All models appear to be in excellent agree- 
ment at least to the 2<x level. A minor difference is the absence 
of lcr PMS solutions in the no- spot case for BCAH98, contrary 
to the CESAM and YY results. 

Additional constraints on the stellar age may be derived 
from CoRoT-2 being a rapid rotator. According to Mamajek & 
Hillenbrand (2008), the 4.5 day rotation period with B - V = 
0.854 (Lanza et al. 2009) implies an age typical of that of the 
Pleiades, i.e., ~ 130 Ma. Given the absence of stars with similar 
B-V and periods shorter than 8 days in the Hyades (~ 625 Ma), 
we estimate that CoRoT-2 is less than 500 Ma old, which thus re- 
stricts the ensemble of solutions from Fig. 5 quite significantly. 

We now focus on the young (< 500 Ma) solutions and com- 
pare CESAM with BCAH98 in the Af*, age) parameter 
space. (We do not show the comparisons with YY, since it 
closely ressembles CESAM). Figure 6 shows that in the no-spot 
case, the lcr solutions are limited to the PMS phase for CESAM, 
and there are no solutions when using the BCAH98 models. At 2 
and 3<r, the solutions span the entire age range and both models 
yield very similar results. We note that MS solutions imply stel- 
lar masses slightly above that of the Sun, whereas PMS solutions 
are distributed between ~ 0.9 and 1.0 M . 

Figure 7 compares the solutions in the (/?*, M*) space. There 
is a clear positive correlation between the two quantities. For 
ages older than 50 Ma, the solutions are confined to high mass 
values and there is a very good agreement between CESAM and 
BCAH98. At young ages however, as seen in previous diagrams, 
the CESAM and BCAH98 solutions differ. The influence of the 
presence of spots is relatively small. 

The results in terms of the stellar mass and radius are sum- 
marized in Table 2 for the different age classes. When no solu- 
tions were found within lcr, the 2cr solutions are indicated. 
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Fig. 6. Constraints obtained for the age as a function of mass 
of CoRoT-2. The left panels correspond to results obtained by 
neglecting the effect of spots. The right panels assume an ad- 
ditional uncertainty in the derived because of the presence 
of up to 20% of spots. The upper panels are calculated from 
CESAM evolution tracks. The lower panels are calculated from 
BCAH98 evolution tracks. Colors have the same meaning as in 
Fig. 5. 
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Fig. 7. Constraints obtained on the mass and radius of CoRoT- 
2. The panels and colors are as in Fig. 6, except that the solu- 
tions are separated between those at young (0 - 50 Ma) and old 
(50 - 500 Ma) ages. Solutions obtained for ages above 500 Ma 
are indicated in grey. 

2.5. Constraints on the planet's physical parameters 

Physical parameters for the planet are derived from the solu- 
tion for the star using the orbital period P or b, the radii ratio 
k = R v /R+, and the semi-amplitude star velocity K (see Table 1) 
using the following equations (e.g. Sozzetti et al. 2007; Beatty 
et al. 2007) 



kR*, 



\2nGl 



smi 



(4) 
(5) 



where G is the gravitational constant, e the orbital eccentricity, 
and i the orbital inclination projected along the line of sight. 

The results in terms of the inferred planet mass and radius are 
shown in Fig. 8 and listed in Table 3. Our inferred planet masses 
are slightly higher than obtained by Alonso et al. (2008) but in 
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Fig. 5. Constraints obtained for the age and radius of CoRoT-2 with different assumptions. The left panels correspond to results 
obtained by neglecting the effect of spots. The right panels assume an additional uncertainty in the derived due to a 0% to 
20% fraction of spots. From top to bottom, the panels are: a) Results obtained with the full CESAM calibrated evolution grid, b) 
Results obtained with CESAM with a mixing length parameter a = O.85a . c) Same as previously but with a = 1.15^ . d) Results 
obtained with the calibrated CESAM evolution grid but a constraint on the stellar density obtained from Alonso et al. (2008) instead 
of Gillon et al. (2010). e) Results obtained with the YY tracks (Demarque et al. 2004). f) Results obtained with the BCAH98 tracks 
(Baraffe et al. 1998). The colored area corresponds to constraints derived from stellar evolution models matching the stellar density 
and effective temperature within a certain number of standard deviations: less than Icr (red), 2<x (blue), or 3<x (yellow). 
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Table 2. Derived mass and radius of the star CoRoT-2 for two different inferred ages. 

Model 



no spot 
MJM Q RJR Q 



0-20% spots 
MJM Q RJR Q 



0-50 Ma 

CESAM 0.97 + 0.04 0.91 +0.02 0.95 + 0.07 0.90 + 0.03 

BCAH98 1 .04 + 0.06 (2cr) 0.93 + 0.03 (2cr) 0.96 + 0.08 0.90 + 0.03 

50-500 Ma 

CESAM 1.03 + 0.01 (2cr) 0.92 + 0.01 (2cr) 1.04 + 0.01 0.93 + 0.01 

BCAH98 1 .03 + 0.02 (2cr) 0.92 + 0.02 (2cr) 1 .04 + 0.02 0.93 + 0.02 

Table 3. Derived mass and radius of the planet CoRoT-2b for two different inferred ages. 



Model 



no spot 0-20% spots 

Mp/M Jup R P /Rjup Mp/M Jup gp/gjup 

0-50 Ma 

CESAM 3.50 + 0.20 1.52 + 0.04 3.45 + 0.27 1.50 + 0.06 

BCAH98 3.68 + 0.25 (2cr) 1.54 + 0.06 (2cr) 3.48 + 0.30 1.50 + 0.06 

50-500 Ma 

CESAM 3.65 + 0.13 (2cr) 1.54 + 0.03 (2cr) 3.67 + 0.13 1.55 + 0.03 

BCAH98 3.65 + 0.16 (2cr) 1.54 + 0.04 (2cr) 3.67 + 0.16 1.55 + 0.04 
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Fig. 8. Constraints obtained on the mass and radius of CoRoT- 
2b. The panels and colors are as in Fig. 7. Rj up = 71492 km 



good agreement with the Gillon et al. (2010) results. The corre- 
sponding planetary sizes are larger than in both studies because 
we account for the effect of spot occultations by the transiting 
planet. The new classes of solutions for the star on the PMS, at 
ages of between 30 and 40 Ma allow however for slightly smaller 
R v values than the standard MS solutions. 



3. Planetary evolution models 

3.1. Modeling procedure 

Our planet evolution models are calculated using CEPAM, a 
code originally derived from CESAM but accounting for the 
physics that is specific to planets (Guillot & Morel 1995). We 
adopt the same two hypotheses as models of other transiting ex- 
oplanets (e.g. Guillot 2008), namely: (i) We use the equation of 
state for hydrogen and helium from Saumon et al. (1995) and a 
slightly larger value of the helium mass-mixing ratio (Y = 0.30) 
to account for the presence of heavy elements; and (ii) interior 
Rosseland opacities are calculated from Allard et al. (2001). The 



outer boundary condition is slightly modified however to allow 
for the determination of the impact of atmospheric properties 
on the planetary evolution. Following Guillot (2010) (see also 
Hansen 2008), we use a T(r) model for the globally-averaged 
temperature field in the atmosphere 



T 4 = 



3T 4 



4 
2 

Yy 



3T e 



eq 



2 

r 



l + 



(6) 



where r is the optical depth at thermal wavelengths, E 2 is the ex- 
ponential integral function of order 2, 4nRpCrTe q is the stellar lu- 
minosity received by the planet, AnR^crT^ is the planet's intrin- 
sic luminosity and y = /c v Mh is the greenhouse factor equal to 
the ratio of the mean visible opacity to the mean thermal opacity. 
In addition, we assume that P = (K t ^/g)r, and link atmospheric 
to interior models at the 10 bar pressure level. 

The values of the coefficients k± and K y are parameter- 
ized from detailed radiative transfer calculations, as described 
in Guillot (2010). These, and hence the atmospheric models in 
general are very uncertain due to the weak constraints on the 
chemical composition of these atmospheres, the unknown cloud- 
coverage, the difficulty in estimating how atmospheric dynam- 
ics transports heat and chemical elements... etc. In any case, the 
values that most closely reproduce the models of Fortney et al. 
(2008) in similar irradiation conditions are k± = 10~ 2 gcm~ 2 
and k y = 4 x 10" 3 gem" 2 . These values yield evolution models 
that closely match our previous calculations (e.g. Guillot 2008). 
However, to match the observational constraints from measured 
brightness temperatures we adopt as a baseline scenario with an 
increased thermal opacity = 1.5 X 10" 2 g cm -2 (see Sect. 3.4). 

The transit radius corresponding to the level for which the 
chord optical depth is equal to 2/3 is calculated as in Guillot 
(2010) using our assumed visible opacity coefficient k y . Because 
of the large size of CoRoT-2b and its relatively significant mass, 
we choose to ignore the possibility that a central dense core 
is present and only consider a fully gaseous planet of approx- 
imately solar composition. An increase in the mean molecu- 
lar weight and/or presence of a core would generally lead to a 
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smaller planet size and thus increase the problem in reproducing 
the observations. 

When considering models that include explicit tidal dissipa- 
tion, we also solve the combined dynamical and structural evolu- 
tion of the star/planet system including tides as described in the 
Appendix. We follow changes in the structural planetary param- 
eters (radius R, pressure P, temperature T, intrinsic luminosity 
L) as well as in the dynamical parameters of the system (semi- 
major axis a, orbital eccentricity e, stellar spin Qi, planetary spin 
Q2). At each evolution timestep, a new value of the tidal heating 
rate H 2 is calculated by solving implicitly the equations for the 
dynamical evolution of the system. This heating rate is used for 
the subsequent structural calculations, and assumed to be dissi- 
pated at the center of the planet, so that L(r = 0) = H 2 . (For a 
discussion of how the depth of the dissipation affects the planet's 
structure and evolution, see Guillot & Showman (2002).) The or- 
bital evolution also modifies the atmospheric boundary condition 
by altering the irradiation flux. The equilibrium temperature T eq 
is hence recalculated at each timestep. Our approach to the tidal 
dissipation calculation is thus similar to that chosen in other cal- 
culations (Ibgui & Burrows 2009; Miller et al. 2009), but it is 
based on the equations derived by Barker & Ogilvie (2009) in- 
stead of those of Jackson et al. (2008). The main differences are 
that the relations are valid for higher values of the eccentricity 
(see also Leconte et al. 2010) and include the secular evolution 
of stellar and planetary spins. 

3.2. Standard evolution models 

We compare in Fig. 9 observational constraints on age and plan- 
etary size to standard evolution models with slightly different as- 
sumptions about mass, initial planetary radius, and helium con- 
tent. The standard models are defined by the planet only con- 
tracting from an initial radius R{ n { as a result of the loss of its 
internal entropy: the only reservoir of energy is the initial grav- 
itational energy J Gm/r dm. The atmospheric boundary condi- 
tion corresponds to our baseline case. 

Our fiducial model has a mass of M = 3.5 Mj, an initial 
radius Ri n i = 2Rj, and an equivalent helium mass mixing ra- 
tio Y = 0.30. It falls short of reproducing the inferred radius by 
~ 20% or more, except at young ages: for pre-main-sequence so- 
lutions at 30 - 40 Ma, the discrepancy is reduced to about 10%. 
Within the error bars, models with different masses have almost 
identical evolutions and therefore planetary mass is not a signif- 
icant parameter in the problem (this is because CoRoT-2b lies 
in the particular regime of mass for which the radius is almost 
independent of mass, which for isolated objects corresponds to 
a maximum in the mass-radius relation -see Guillot (2005)). In 
a similar way, changing the initial radius affects the evolution 
only in the first million years: the initial conditions are rapidly 
forgotten. One may wonder whether a different composition, in 
particular a different helium abundance, would have a stronger 
effect on evolution? As shown in Fig. 9, models with Y = 0.25 
indeed yield a ~ 4% larger size than for Y = 0.30, but this again 
falls short in explaining the large size of the planet. 

These calculations confirm that CoRoT-2b is an anomalously 
large planet, a result already obtained by Alonso et al. (2008), 
Leconte et al. (2009), and Miller et al. (2009). However it also 
demonstrates the fact that the planet's young age is likely to be a 
crucial factor in explaining its size, both because of the possibil- 
ity that the planet is initially quite large, and because our stellar 
evolution models yield solutions at 30-40 Ma that are closer to 
the theoretical evolution tracks than at any later times. 




1.0 1 ^ i 

0.001 0.010 0.100 1.000 10.000 
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Fig. 9. Contraction of CoRoT-2b compared to its measured ra- 
dius and inferred age. Our fiducial model has an initial radius 
Ri n i = 2Rj, mass M p = 3.5 Mj, and equivalent helium mass- 
mixing ratio Y = 0.30. The evolution for models with different 
initial radii Ri n i (3 to 6Rj), different masses (3.3 to 3.7 Mj), or a 
different value of Y (0.25) are shown as labeled. 



3.3. CoRoT-2b among its peers 

It is instructive to compare CoRoT-2b to an ensemble of other 
transiting giant planets. Among these, CoRoT-2b may not be the 
largest (it is smaller than CoRoT-lb, HAT-P-8b, TrES-4b, and 
the present record-holder WASP- 12b), but it remains the most 
difficult to reconcile with present-day models. This is most easily 
shown by ranking the planets in terms of their radius anomaly, 
i.e. the difference between its measured radius and the one pre- 
dicted by models of a pure solar-composition planet of the same 
mass and age (Guillot et al. 2006). As shown by Fig. 10, when 
using standard models, CoRoT-2b has a positive, large radius 
anomaly of 20000 km, but still smaller than that of HAT-P-8b, 
TrES-4b, and WASP- 12b. However, for these last three planets, 
their large radius can be explained within the error bars by an 
additional heat source equivalent to 1% of the incoming stellar 
luminosity deposited at the planet's center (see Guillot 2005 and 
references therein). As shown in the right panel of Fig. 10, this 
is not true for CoRoT-2b: because of its large mass, the planet 
tends to contract rapidly and therefore requires special condi- 
tions to explain its large size. 

3.4. CoRoT-2b's atmosphere 

We now consider how models of the planet's atmosphere affect 
its evolution. Remarkably, secondary transits of CoRoT-2b were 
detected in the optical from C0R0T lightcurves (Alonso et al. 
2009; Snellen et al. 2010) and in the infrared from Spitzer IRAC 
observations (Gillon et al. 2010) as well as ground-based mea- 
surements (Alonso et al. 2010a). These directly probe the plane- 
tary atmosphere and are thus key constraints of the outer bound- 
ary conditions used in the evolution modeling. 

The fluxes inferred from these measurements are equiva- 
lent to brightness temperatures of 1325 ± 180K at 8 /mi and 
1805 ± 70 K at 4.5^m (Gillon et al. 2010; Snellen et al. 2010). 
Additional ground-based measurements yield 7b = 1890+^ K 
in the band (~ 2.1 jim) and an upper limit of 2250 K in the H 
band (~ 1.6/zm) (Alonso et al. 2010a). In the optical, indepen- 
dent studies yield brightness temperatures that are very similar 
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Fig. 10. Radius anomaly (difference between observed and mod- 
eled radius -see text-) as a function of planetary mass (in log) 
for a selection of known transiting giant planets. Left panel: 
Standard models (including stellar irradiation but no extra heat 
source) are used for the comparison. Right panel: Models as- 
sume an extra heat source at the planet's center equivalent to 
1% of the incoming stellar heat flux (see Guillot & Showman 
2002). CoRoT-2b is labeled "27". It is the most anomalously 
large planet on the right panel. 



within error bars, i.e. 2120+*°, K (Alonso et al. 2009, 2010b) 
and 2170 ± 50 K (Snellen et al. 2010). A crucial consequence 
that can be derived is that independently of hypotheses about 
the atmospheric composition and opacity sources, the day-side 
atmosphere of CoRoT-2b is characterized by physical tempera- 
tures that are as low as 1325 ± 180^, and at least as high as 
1805 ± 70 K. The case of the optical brightness measurements 
are more complex to interpret because they may arise from ei- 
ther direct emission or a reflection of the incoming stellar flux 
(Alonso et al. 2009; Snellen et al. 2010). In the limit of a geo- 
metric albedo of 0.2, the flux would be entirely due to direct re- 
flection and thus provide no information about the atmospheric 
temperature profile. In contrast, a zero albedo would imply that 
the flux in the optical is thermal emission from the planet, and 
that the corresponding temperatures are high. 

These temperature constraints are compared in Fig. 11 to 
temperature profiles calculated for CoRoT-2b in the framework 
of our semi-gray analytical model assuming a full redistribution 
of the incoming stellar flux (see Eq. (6) and Guillot 2010). The 
value of the intrinsic flux 7^ = 1000 K was chosen to match 
that of models with a size ~ 1.57?j up , as observed. Using values 
of the thermal and visible opacities calibrated to detailed atmo- 
spheric models (Fortney et al. 2008) (labeled /c* h , k* in the figure), 
we derive a temperature profile that is difficult to reconcile with 
the Spitzer and C0R0T brightness temperatures: the temperature 
range spanned by the atmosphere at mean optical levels smaller 
than unity is small: from about 1400 to 1600 K. This implies that 
to explain the observations with this model, one needs to invoke 
(i) an albedo of ~ 0.2, (ii) a high-opacity at 8 yum, and (iii) a very- 
low opacity (lower than at visible wavelengths) at 4.5 fim. This 
combination of factors appears to be unlikely. Furthermore, one 
should consider that the low 8 jim temperature has to be emit- 
ted from high levels in the atmosphere. Because the radiative 
timescale is shorter at high altitudes (Showman & Guillot 2002; 
Iro et al. 2005), its temperature may have to be closer to the day- 
side than global average, which would make the problem even 
worse. 
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Fig. 11. Possible atmospheric pressure-temperature profiles for 
CoRoT-2b compared to observational constraints. The bright- 
ness temperatures obtained by Spitzer IRAC at 4.5 and 8//m 
(Gillon et al. 2010) and by C0R0T in the optical (Alonso et al. 
2010b; Snellen et al. 2010) are indicated by vertical grey bands. 
In the case of the optical brightness, the corresponding physical 
atmospheric temperatures depend on the atmospheric geomet- 
ric albedo (A = to 0.15), and are indicated by squares and 
with error bars. Temperature profiles are calculated on the ba- 
sis of a semi-grey analytical model (Guillot 2010), with fiducial 
values of the thermal and visible opacities K* h - 0.01 cm 2 g" 1 , 
< = 0.004 cm 2 g _1 , and T mi = 1000 K. The different lines cor- 
respond to different values of these coefficients (as labeled). The 
values of the assumed Bond albedo are A = (purple, orange, 
red curves), A = 0.1 (yellow curve) and A = 0.15 (blue curves). 
The levels for which thermal and visible optical depths equal 2/3 
are indicated by diamonds and triangles, respectively. 



For these reasons, we consider alternative models. We re- 
strict ourselves to models without a temperature inversion, both 
because this is not predicted by dedicated radiative transfer mod- 
els (Gillon et al. 2010; Snellen et al. 2010), and it would make 
the planetary radius problem more severe. Figure 1 1 presents a 
variety of alternative profiles. The temperature range in the low- 
optical thickness part of the atmosphere is directly related to the 
factor y = k w /k±, with a low T value implying a strong green- 
house effect. We find that models that can more readily repro- 
duce both the infrared and visible brightness temperatures are 
those with y « 0.2 to 0.4. These models are also consistent with 
the secondary transit of CoRoT-2b being detected in CoRoT's 
red channel and not in the blue and green channels (Snellen et al. 
2010), implying that it indeed originates mostly from direct plan- 
etary emission rather than stellar light reflection. They are also 
consistent with their smaller implied albedos, in line with the low 
limit A g = 0.038 ± 0.045 obtained for HD209458b (Rowe et al. 
2008), a planet with a similar equilibrium temperature. Larger 
greenhouse effects (i.e. the red curve in fig. 11) are not consis- 
tent with the models because they would tend to yield brightness 
temperatures that are larger than inferred. 

Compared to traditional atmospheric models, this larger 
greenhouse effect may be achieved by several means: the pres- 
ence of clouds, of photochemical products, or generally of minor 
species that are efficient infrared absorbers but are transparent at 
visible wavelengths. Complications arising from chromatic and 
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Fig. 12. Transit radius of CoRoT-2b as a function of age. The 
colored areas indicate constraints obtained from photometric, 
spectroscopic, radial- velocity data and stellar evolution mod- 
els. The colors are a function of the distance to the effective 
temperature- stellar density constraint ellipse in standard error 
units: within l<x (red), 2<x (yellow) and 3<x (blue). The star's 
activity implies an age smaller than ~0.5 Ga, which indicated 
by solutions in shades of grey. The main constraint is obtained 
for a lightcurve analysis that includes the presence of spots. An 
analysis that does not account for spots yields a ~ 3% smaller 
radius, indicated by a dotted area (see text). The plain curves 
are obtained from standard evolution models using atmospheric 
boundary conditions that are parameterized by mean infrared k± 
and visible /c v opacity coefficients, as in Fig. 11. The dashed and 
dotted curves are models from Leconte et al. (2009) and Gillon 
et al. (2010), respectively. 



dynamical effects should of course be taken into account, but are 
should not change qualitatively the conclusions of this work. 

3.5. Impact of atmospheric models on the planetary evolution 

Planetary evolution models calculated with different atmo- 
spheric boundary conditions are compared in Fig. 12 to standard 
evolution models for this planet available in the literature and 
to our inferred age/size constraints for CoRoT-2b. We first note 
that our baseline models (with opacities K* h and k*) are a close 
match to the models of Gillon et al. (2010) after about 20 Ma, 
and that they reproduce relatively well the models of Leconte 
et al. (2009), within 5%, over the whole age range. In this last 
case, the differences in behavior may be attributed to our fidu- 
cial model having slightly lower infrared opacities and a smaller 
y value. Because the internal structure part of the calculation 
should be very similar, this highlights that even with similar hy- 
potheses (a solar composition cloud-free atmosphere), detailed 
radiative transfer atmospheric models yield predictions that dif- 
fer in a relatively significant way (see also Guillot 2010). 

As previously, most models are unable to reproduce the con- 
straints. However, two of them do intersect the constraint area at 
young ages between 30 and 40 Ma. These are models for which 
the thermal atmospheric opacity has been increased by a fac- 
tor 3.5 to 4, relatively independently of visible opacities. The 
consequence of the larger k± is to reduce the amount of initial 
heat lost by the planet at young ages. The intersection is small, 
but could become larger if we considered that non-grey mod- 
els may enable the deep atmospheric temperatures to increase 



Fig. 13. Contraction of CoRoT-2b relative to its measured radius 
and inferred age (see fig. 12). The evolution models are labeled 
as follows: standard (black): irradiated, solar-composition planet 
with no extra heat source; 1% K.E. (blue): 1% of the incoming 
stellar flux is assumed to be dissipated at the planet's center; 
opacitiesx30 (light blue): opacities have been artificially multi- 
plied by 30 compared to the standard model; 10 29 (red): model 
in which 10 29 erg s" 1 is deposited at the planet's center. 



while the visible brightness temperature remains constant, and 
also that the planet probably formed at least a few million years 
after the star. As shown in Fig. 12, for an even stronger green- 
house effect (a smaller y), the present size of the planet may be 
reached at even older ages. However, this extreme model is un- 
likely because it disagrees with the inferred visible brightness 
temperature (see § 3.4). 

The present size of CoRoT-2b can thus be explained by the 
combination of a young age of between 30 and 40 Ma, and ad- 
ditional opacity sources (gases/clouds) in the atmosphere. We 
note that Burrows et al. (2007) also proposed an increase in at- 
mospheric opacities to explain the large sizes of exoplanets. Our 
solution is similar, but we emphasize that this increase should 
concern more particularly the opacities at infrared wavelengths. 



3.6. Alternative recipes 

Alternative models invoked to explain the large sizes of other ex- 
oplanets are unlikely to work. As shown in Fig. 13, the kinetic- 
energy dissipation model proposed by (Guillot & Showman 
2002) and used successfully for most known transiting planets 
(Guillot et al. 2006; Guillot 2008) fails for CoRoT-2. This is also 
the case of a considerable -30 fold- increase in interior opac- 
ities that would also explain the sizes of most transiting plan- 
ets (Guillot 2008). In fact, as already noted (Alonso et al. 2008; 
Gillon et al. 2010), the energy dissipation deep inside the planet 
required to explain the present-day radius is enormous, on the 
order of 10 29 ergs -1 . This is about 30000 times the present in- 
trinsic luminosity of Jupiter. It is also about l/4th of the power 
that the planet receives from its parent star. 

After that provided by stellar radiation, the most important 
potential source of energy is that taken from the planetary or- 
bit. When moving CoRoT-2b from infinity to its present orbit, 
AE = GM*M v /2a « 10 45 erg have to be dissipated. If this energy 
dissipation were to occur entirely in the planet, the maximum 
amount of time one would be able to maintain a 10 29 erg s" 1 dis- 
sipation rate is ~ 300 Ma. 
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3. 7. The effect of tides 

As originally proposed by Bodenheimer et al. (2001) and Gu 
et al. (2003) and later studied by many authors (e.g. Jackson et al. 
2008; Ibgui & Burrows 2009; Miller et al. 2009), stellar tides 
provide a way to transfer gravitational energy from the plane- 
tary orbit into the planet and either slow its contraction, or even 
produce a size inflation. Models coupling the equations govern- 
ing the dynamical evolution of the star+planet system with the 
physical planetary evolution rely however on a crucial assump- 
tion: that dissipation occurs at a sufficient depth in the planet 
interior, i.e. roughly within the planet's convective zone, deeper 
than a few 100 bars or so (see Guillot & Showman 2002, for a 
discussion). The mechanisms responsible for the dissipation are 
yet unknown, and may occur either high up in the atmosphere 
(Lubow et al. 1997) or throughout the planetary interior (Ogilvie 
& Lin 2004). 

Following Gillon et al. (2010), we present models of the dy- 
namical and physical evolution of the CoRoT-2 system caused 
by the action of stellar and planetary tides. We maximize the 
efficiency of the heat dissipation by assuming that it is entirely 
deposited at the center of the planet. We use the dynamical evo- 
lution equations derived by Barker & Ogilvie (2009) and include 
high order terms in eccentricity and equations for the evolution 
of the stellar and planetary spin (see Appendix). On the basis of 
the calculations by Jackson et al. (2008), we explore values of 
the tidal factor g p between 10 5 and 10 6 , and of g* of 10 5 and 
higher. 

We analyze in Fig. 14 how the tidal heating rates and the 
orbital timescales depend on the eccentricity of the system us- 
ing all other known parameters of the system. We first note that 
for values of the eccentricity e > 0.3, these become extremely 
stiff functions of e. This implies that any initially high eccen- 
tricity value causes a rapid evolution of the system which is 
hardly predicted by models developed only to second order in 
eccentricity (Jackson et al. 2008; Ibgui & Burrows 2009; Miller 
et al. 2009; Gillon et al. 2010). We find that a 10% asynchronous 
planet would dissipate the required luminosity, but the corre- 
sponding synchronization timescale is extremely short, about 
10, 000 years. 

At low eccentricities, an inward migration with a timescale 
~ 1 Ga(g*/10 6 ) results from tides raised by the planet onto the 
star.At high eccentricities and for our choice of Q factors, tides 
raised by the star onto the planet begin to dominate and cause a 
decrease in the semi-major axis that is concomitant to the circu- 
larization of stellar tides on the planet. The orbit circularization 
is mostly caused by the planet, unless Q p > 10<2*. While the 
planet is synchronized efficiently, the star is found to be spun up 
by the planet relatively slowly ~ 1 Ga for eccentricities e < 0.2. 

Given its inferred eccentricity, the present size of CoRoT- 
2b may be explained by tides only within two scenarios: (i) By 
a very low Q v value and a forced eccentricity due to the pres- 
ence of another planet, (ii) By an initial stage of high-dissipation 
followed by a rapid circularization and contraction. The former 
case is unlikely (see also Gillon et al. 2010). The last possibility 
requires that the circularization proceeds faster than the planet's 
contraction. 

In Fig. 15 we explore the constraints that can be derived on 
Q p . Using the equations in the Appendix, we calculate the min- 
imum eccentricity required for tides to dissipate 10 29 ergs _1 in 
the planet (top panel). We then calculate the time required for the 
eccentricity to decline from this value to the observed one (we 
assumed e = 0.02). This time is compared to the time required 
to contract from 2 to 1.5 Rj based on our different atmospheric 
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Fig. 14. Heating rates and dynamical timescales for the orbital 
evolution of the CoRoT-2 system as a function of the orbital ec- 
centricity. The star and planet are assumed to have their present 
mass and orbital period, and tidal dissipation factors g* = 10 6 , 
<2p = 10 5 . The blue curves refer to the star, the red curves to the 
planet. The planet is assumed to be either synchronously rotating 
(spin period=orbital period) (plain lines), or to rotate 10% faster 
than synchronous rotation (dashed lines). From top to bottom, 
the panels show various important quantities: (a) Heating rates 
due to tidal dissipation in the planet and in the star, (b) Migration 
timescale due to the star and the planet. Migration is inward ex- 
cept in the asynchronous case due to the planet's spin down for 
eccentricities smaller than ~ 0.12. (c) Spin-up timescale: The 
star is spun up by the planet in all cases. The planet is gen- 
erally spun up, except in the asynchronous case (see text), (d) 
Circularisation timescale: The evolution towards a circular orbit 
is due to both stellar and planetary tides. 



models (bottom panels). We thus derive an upper limit to the 
planetary tidal factor Q v ~ 10 6 that is similar to the 10 5 5 quoted 
by Gillon et al. (2010). This estimate does not account for the 
effect of migration but we propose that it should be relatively 
realistic. 

We show in Fig. 16 a few example results of the full dy- 
namical calculation. One problem we faced was the existence of 
a runaway inflation of the planet especially for high values of 
the eccentricity. In order to maximize the ensemble of solutions 
(and ease the calculation), we set an arbitrary saturation value 
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Fig. 15. Minimum eccentricity and circularization timescale as 
a function of tidal dissipation factor Q' p . Top panel: Eccentricity 
required to reach a dissipation rate H = 10 29 erg s" 1 in the planet. 
Bottom panel: Time to evolve from the minimum eccentricity to 
the observed e ~ 0.02 for two values of the stellar tidal factor, 
oo (dashed) and 10 6 (plain). The horizontal lines show the times 
required for CoRoT-2b to contract from 2 to 1.5 Rj, with the at- 
mospheric boundary conditions as in Fig. 1 1 . 



at H = 2 x 10 29 erg s" 1 and increase the planetary Q p factor ac- 
cordingly. Using this recipe, we are able to reproduce approxi- 
mately the solutions found by Gillon et al. (2010), but find that 
the contraction is too fast compared to our age constraints. We 
were unable in general to find solutions when starting from a 
fixed initial eccentricity at time t = 0. We can obtain transient 
solutions instead when including inward migration or a late in- 
crease in the eccentricity. A plausible scenario could thus be that 
CoRoT-2b had a recent (< 20 Ma) planetary encounter leaving 
it into an eccentric orbit (e.g. Juric & Tremaine 2008; Ford & 
Rasio 2008). Another similar possibility is that the planet had 
a Kozai interaction with another distant body (e.g. Fabrycky & 
Tremaine 2007), its orbit was changed to one of high eccentric- 
ity, and circularization began less than 20 Ma ago. 



3.8. A recent planetary impact? 

Another alternative is of course that CoRoT-2b experienced a 
planetary impact in the past ~ 20 Ma. Juric & Tremaine (2008) 
find that up to about 20% of planets may experience mergers 
in systems with multiple embryos of the same mass. This prob- 
ability is extremely dependent on initial conditions, but illus- 
trates that giant impacts may not be an extremely rare possibil- 
ity. We illustrate the possible mass of such an impactor with an 



extremely simple model. The energy required to inflate the pre- 
cursor of CoRoT-2b is 

GM 2 



(7) 



o 



where Mo is the mass of the precursor, Ro its size, and is the ra- 
dius factor increase required to explain CoRoT-2b's present size. 
We assume that the collision takes place at the escape velocity 
of the precursor, and that half of the kinetic energy is transferred 
into increasing the internal energy of the final planet, i.e., that 

1 GMqMi 

Azi collision ~ 



Ro 

where M\ is the mass of the impactor. Using M p 
as the present mass of CoRoT-2b, and AE co ui S ion ' 
obtain 



Mi 



(8) 
- M (} + M\ 

AEg rav We 

(9) 



This relation is independent of the initial planetary radius. This 
is because both the impact energy and the inflation energy have 
the same Ro dependency. If we assume ft « 1/4, we find 
Mi ~ 1/2 M p ~ 1.75 Mj. If the system is really young, the need 
for a large impactor decreases: e.g. with j3 « 0.05, we obtain 
Mi ~ 1/6 M p ~ 0.6 Mj. These values of course represent ex- 
tremely simplified examples but they show that unless planets 
can be accelerated to much higher velocities (this would prob- 
ably require a fourth, massive high density object, i.e. a brown 
dwarf), only collisions between giant planets would provide the 
necessary energy to significantly inflate CoRoT-2b. 

4. Conclusion 

We have combined stellar and planetary evolution models to help 
us develop a consistent scenario to understand the formation and 
evolution of the CoRoT-2 system. Although stellar spots are im- 
portant in the lightcurve analysis and do have an important im- 
pact when deriving the planetary radius, we have demonstrated 
that their consequences on the derived stellar properties are mod- 
est and can be modeled as an additional uncertainty in the de- 
rived effective temperature of the star. 

We have presented evidence for the youth of the CoRoT- 
2 system. The rapid ~ 4.5 day spin of the star is of course a 
strong indication of a young age, but we have also shown that 
a very promising ensemble of solutions occurs on the pre-main- 
sequence phase of the parent star, yielding ages of between 30 
and 40 Ma. By combining constraints obtained on the atmo- 
sphere from both brightness temperature measurements and stel- 
lar and planetary evolution models, we found that three scenarios 
can explain the present large size of planet CoRoT-2b. All these 
imply a recent event that took place less than 40 Ma ago. 

1 . Its atmosphere is about 3 to 4 times more opaque in the in- 
frared than usually thought, and the system is indeed 30 to 
40 Ma old. 

2. The CoRoT-2 system involved multiple planets and less than 
20 Ma ago, close encounters between the planets brought 
two into collision. The CoRoT-2b impactor would need to 
have had a relatively significant fraction of the total plane- 
tary mass. 

3. The CoRoT-2 system involved multiple planets and less than 
20 Ma ago, close encounters between the planets left CoRoT- 
2b on an eccentric orbit. This orbit has been almost circular- 
ized, but remains in the transient heating phase. The parame- 
ter space required in terms of initial eccentricities and orbital 
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Fig. 16. Evolution of the CoRoT-2 system in the presence of tides, as a function of age expressed in billion years. From top to 
bottom, the five panels are: (a) Planet's transit radius; (b) Planet's orbital eccentricity; (c) Planet's orbital period in days; (d) Planet's 
spin period in days; (e) Star's spin period in days. The observational constraints are indicated by colored areas (see fig. 12). For 
visibility purposes, the constraint on the orbital period has been increased by a factor 10 5 . Three models are shown: Black lines: 
A model with Q v = 3 x 10 5 , g* = 10 5 , e ini = 0.4, a ini = 0.04 AU, corresponding to a preferred model of Gillon et al. (2010). 
Purple lines: A model including migration for the first 10 Ma, and Q v - 10 6 , = 10 6 , e^ = 0.4, = 0.13 AU. Blue lines: A 
close-encounter model, for which the eccentricity is suddenly raised to e ~ 0.77 after 300 Ma, presumably due to close encounters 
with another planet. This assumes Q' p = 10 5 , g* = 10 6 > a mi = 0.1 AU. The thick parts of the lines corresponds to epochs for which 
the eccentricity and orbital periods are in agreement with the observations (see text). 
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distances for this to have occurred is small, and we found that 
it generally requires a saturation of the tidal heating to avoid 
the complete loss of the planet. 

The detection of abundant lithium in CoRoT-2 confirms that 
the star is young: the comparison with data obtained for open 
clusters indicate that its age should be between 30 and 320 Ma 
(Gillon et al. 2010). However, this may be an overestimation be- 
cause lithium may be more easily destroyed in the presence of a 
massive protoplanetary disk (see Bouvier 2008). Several kinds of 
observations of CoRoT-2 would shed light on its nature. The de- 
tection of an infrared excess or a debris disk would be an indica- 
tion that the system is young even though there is a wide spread 
in disk ages (e.g. Hillenbrand et al. 2008). Searching for distant 
companions would also help especially given the possibility that 
the planet's eccentricity may have been pumped up by Kozai in- 
teractions with a third body before being efficiently damped by 
tides. Last but not least, the determination of a complete infrared 
lightcurve including both the primary and the secondary tran- 
sit would be extremely valuable to constrain the planet's atmo- 
spheric properties, in particular its day-night heat redistribution 
efficiency. It would also help us to determine the planet- to- star 
radius ratio in way that is less affected by systematic errors due 
to stellar activity. With similar observations in the visible, one 
would be able to test our assumption that the signal detected in 
the CoRoT lightcurves is due to thermal emission from deep lev- 
els in the planet. Since the radiative timescales rapidly increase 
with depth (Iro et al. 2005), we would expect that phase varia- 
tions in the visible are smaller than in the infrared. 
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Appendix 

We present the time-averaged equations for the dynamical evo- 
lution of the star-planet system used in this work. These were 
taken from Barker & Ogilvie (2009) and applied to the case of 
a planar system (i.e. neglecting possible inclinations of the star 
and planet), accounting for external forcing of the semi-major 
axis, eccentricity, and the conservation of angular momentum 
during the planet's contraction. 

Given a semi-major axis a and an eccentricy e, the orbital 
mean motion is n = (G(m\ + m2)/a 3 ) l/2 and the angular momen- 
tum h = na 2 (l - e 2 ) 112 . The equations of secular evolution of h, 
e, the star's spin Q\ and planet's spin Q2 under the effect of tides 



are: 
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are 



In these equations, fi = m\m 2 /(m\ + m 2 ) is the reduced mass 
of the system, I\ 2 are the moments of inertia of the star and 
planet, and tf\ r2 are the tidal friction timescales, defined as 
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where Q[ 2 are tidal dissipation efficiency parameters defined as 
Q\ 2 = 3/(2^i ? 2^Ti ? 2), ^1,2 being the second-order Love numbers 
for the star and planet, and T\^ 2 the assumed constant lag time be- 
tween the quasi-hydrostatic figure and its equilibrium tide value 
(see Barker & Ogilvie 2009; Eggleton et al. 1998, for a discus- 
sion). 

The quantities a, e, Cl\, and Cl 2 correspond to imposed rates 
of change per unit time of the semi-major axis, eccentricity, star 
spin, and planet spin, respectively. They are all set to mimic the 
presence of additional physical processes: a ± accounts for the 
migration imposed by a circumstellar disk on the planet, e 
mimics the mean eccentricity increase due to a Kozai effect, and 
£l\ ± accounts for the evolution of the stellar spin due to mag- 
netic breaking. The rate of change of the planet spin Cl 2 however, 
is self-consistently obtained from angular momentum conserva- 
tion considerations that account for the planet's contraction rate 
R 2 obtained from the evolution calculations 



Q 2 . 

Q 2 = -2—R 



(16) 



This relation assumes that the planet rotates as a solid body, and 
that its moment of inertia remains constant. 

Furthermore, Eqs (10-13) make use of the following func- 
tions of the eccentricity (Barker & Ogilvie 2009; Hut 1981): 
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Last but not least, from energy conservation, one infers that 
the energy dissipated in the star and in the planet, respectively, 
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The heating rate in the planet H 2 may then be applied to the 
planetary evolution calculations. (Note that Eqs. (23) and (24) 
have opposite signs compared to Barker & Ogilvie (2009), so 
that Hi and H 2 are positive). 

For models in which tidal dissipation in the planet is assumed 
to saturate at a fixed value H 2 = H 2yI{ax , we proceed as follows: 
when the dissipation in the planet is found to exceed H 2 ^ max , we 
increase the value of Q' 2 until we reach the maximal allowed 
dissipation. 
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